First principles description of the giant dipole resonance in ^^O 
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We present an ab-initio calculation of the giant dipole resonance in ^^0 based on a nucleon- 
nucleon (NN) interaction from chiral effective field theory that reproduces NN scattering data with 
high accuracy. By merging the Lorentz integral transform and the coupled-cluster methods, we 
extend the previous theoretical limits for break-up observables in light nuclei with mass numbers 
(j4 < 7), and address the collective giant dipole resonance of ^^0. We successfully benchmark the 
new approach against virtually exact results from the hyperspherical harmonics method in *He. 
Our results for reproduce the position and the total strength (bremsstrahlung sum rule) of the 
dipole response very well. When compared to the cross section from photo-absorption experiments 
the theoretical curve exhibits a smeared form of the peak. The tail region between 40 and 100 MeV 
is reproduced within uncertainties. 
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Introduction. — Giant resonances dominate the dipole 
response in a variety of self-bound quantum systems such 
as atomic nuclei [1 , atomic clusters |2], and fuUerenes [3]. 
In photonuclear cross sections, the giant resonance is a 
broad and very pronounced peak at energies around 10- 
20 MeV. After its first observation in Uranium by Bald- 
win and Klaiber in 1947 [3], this prominent structure 
has been observed across the nuclear table. The first 
theoretical interpretations were given by Goldhaber and 
Teller [5], who proposed that the observed peak had to 
be of dipole nature [hence the name giant dipole reso- 
nance (GDR)], and by Steinwedel and Jensen [B]. It was 
described as a nuclear resonance due to the collective mo- 
tion of protons against neutrons, whose width was due 
to the transfer of energy from this vibration into other 
modes of the nuclear motion [7]. Since then the GDR 
has been the object of numerous theoretical studies that 
have tried to account for the GDR centroid and width in 
collective as well as microscopic many-body approaches. 
The microscopic computation of giant dipole resonances 
is a challenging task. In self-consistent mean field theo- 
ries, efficient calculations employ finite amplitude meth- 
ods 0, and the microscopic description of the GDR's 
centroid was only achieved very recently by including the 
quasiparticle-phonon coupling ^ . 

First principles computations of the GDR are only 
available for very light nuclei with mass numbers A < 7. 
One of the main challenges in such computations is that 
the GDR appears in the continuum energy region. Many- 
body approaches not only need to account for compli- 
cated ground-state correlations, but in principle also need 
to solve the many-body scattering problem. In ab-initio 
few-body calculations these challenges were overcome by 
employing the Lorentz integral transform (LIT) [TOl [H] . 
This method reduces the continuum problem to a bound 



state- like problem and thereby overcomes one of the main 
challenges. The nuclear few-body problem is solved with 
the hyperspherical harmonics (HH) expansion [121 US] 
and the no-core shell model [Tl]. Starting from semi- 
realistic NN interactions, this approach computed the 
"resonant shape" of the photonuclear cross section of 
^He HH Hg, 6Li/6He [HI HHj and ^Li [H]. For ^He, 
modern realistic two- and three-body forces were also 
employed in the computations [23 HI] ■ 

The studies of the A = 6 isotopes [ITl [18] showed that 
the halo structure of the rare isotope ^He differs consid- 
erably from the stable ®Li nucleus. In ^He, the dipole 
cross section consists of two well separated peaks, while 
a single resonant shape is observed for ^Li. These re- 
sults do not depend on the employed NN interaction. 
It would be interesting to extend the ah initio studies 
of the nuclear electric dipole response to heavier exotic 
nuclei such as neutron-rich isotopes of Oxygen and Cal- 
cium. The dipole cross section of exotic Oxygen isotopes 
exhibit fragmented low-lying strengths and thereby dif- 
fer from the stable ^^O [22l [23] . As a first step toward 
this goal it is necessary to extend the applicability of the 
ab-initio LIT method from very light nuclei to heavier 
isotopes. In this Letter, we will combine the LIT with 
the coupled-cluster (CC) method [24l[25], and microscop- 
ically compute the GDR for ^^O. Our calculations employ 
NN interactions from chiral effective field theory [551 127] ■ 

Method. — The inclusive response function, also called 
dynamical structure function, to a dipole excitation, is 
defined as 

Siio) = J2 li^olDoim'^Ef ^Eo-Lo). (1) 

Here uj is the excitation energy, Jq is the nuclear ground- 
state spin, while ^pQ and ipf denote the ground state and 
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final state wave functions of the nucleus with energies 
Eq and Ef, respectively. The translationally invariant 
dipole operator is 

D-^P.(r,-R,„0=^('p,-|')r,. (2) 

i i ^ ^ 

Here and Rem are the coordinates of the i-th particle 
and the center-of-mass, respectively, while Pi is the pro- 
ton projection operator and Z is the number of protons. 
The zero-component of the dipole appearing in Eq. ([T]) 
is Dq = D • Gq, where Bq is a unit vector pointing in the 
direction of the photon propagation. 

The symbol in Eq. (jl| indicates both a sum over 
discrete energy levels and an integration over continuum 
eigenstates. For nuclei with mass numbers A > 3 one 
faces the problem that when -0/ is in the continuum, it 
cannot be calculated in an exact way at every excitation 
energy. Thus, a direct calculation of S(uj) is not possible, 
unless one introduces approximations. However, via the 
LIT method the problem can be reformulated in such a 
way that the knowledge of -0/ is not necessary. To this 
end, one introduces the LIT of the response function as 

L(u;o,r) = rck.- ^^^^—^ ■ (3) 

Here wth is the threshold energy and F > 0. By using 
the closure relation one finds 

Liz) = {^o\Dl ~^^ Do\^o) = {m, (4) 

H ~ z* H ~ z 

where we introduced the complex energy z = Eo+uJo+iT. 
The LIT of the response function, L(z), can be computed 
directly by solving the Schrodinger-like equation 

[H - zM ^ DoHo) (5) 

for different values of z. The solution ijj has bound-state- 
like asymptotics and does not pose the difficulty that 
one encounters in the direct calculation of the continuum 
eigenstates ipf. Thus, L{z) can be calculated even for 
A > 3 with any good bound-state techniques, and results 
for the strength function S{u;) are typically obtained from 
a numerical inversion pSVf^U] of the LIT ([s]). 

We are interested in the computation of the dipole re- 
sponse for ^^O. To this end, we merge the LIT with the 
CC method 'MJ |551 [3T] , which is a very efficient bound- 
state technique, applied with success on several medium- 
mass nuclei [321 - 136] . In CC theory the ground state of 
the system is given by |?/'o) — exp(T)|0o)- Here |0o) is a 
Slater determinant and T generates particle-hole (ph) ex- 
citations. One has T = Ti+T2 + . . ., as the sum of a Ip-lh 
operator Ti, a 2p-2h operator T2, ect. In what follows we 
will consider only an expansion up to T2, which is known 
as coupled-cluster with singles-and-doubles (CCSD) exci- 
tations ^25]. In CC theory the LIT of the dipole response 
function is 

L{z) = {Ol\dI ^[ Do\Or) = ii^Rli^L) . (6) 



Here H = exp(— T)i7 exp(r) is the similarity trans- 
formed Hamiltonian, and D^P — exp(— T)Z)q^-' exp(T) is 
the similarity-transformed dipole (dagger) operator. The 
states (Oil and |0i^) are the left and right ground-states 
of the non-hermitian similarity-transformed Hamiltonian 
H, respectively. In analogy to Eq. we compute 
the LIT of the response function by solving a right 
Schrodinger-like equation 

(H-z)IV^fl(z)) =^o|Ofl,), (7) 

and an equivalent left Schrodinger-like equation for 
'0l(^). Here, 10^(2:)) = R{z)\(j)o), and the operator R 
is a linear expansion in particle-hole excitations, 

R{z) ^Ro + J2 ^^^^ + ^j^^^A + ■■■ ■ (8) 

ia ijab 

To be consistent with the CCSD approximation we trun- 
cate R{z) at the 2p-2h excitation level. It is implied 
that sums over i,j (a, &) run over occupied (unoccupied) 
states. An equivalent expansion can be written for the 
left excitation operator. The LIT of Eq. ([6| can be com- 
puted efficiently by employing a generalization of the 
Lanczos algorithm for non-symmetric matrices in the so- 
lution of the linear problem ([7| . Details will be shown in 
a forthcoming publication. In our calculations we employ 
a NN interaction from chiral effective field theory at next- 
to-next-to-next-to leading order [37j, supplemented by a 
point Coulomb force. We omit three-nucleon forces that 
already appear at next-to-next-to leading order |38| due 
to the problems associated with computing three-nucleon 
forces in the large model spaces we employ. Note that 
three-nucleon forces have a smaller impact on the photo- 
absorption cross section of *He [^Hj (few percent in the 
peak), but it is not clear how large their effect is for -'^^O. 

Benchmark for ^He. — For quality assurance, we want 
to compare the CC results with virtually exact results ob- 
tained from the effective interaction hyperspherical har- 
monics (EIHH) method 112^ for ''He. Figure [ija) shows 
L{ujo,T) as a function of ujq at fixed F = 10 MeV, 
computed with the CC method in model spaces with 
nn = 20 MeV and iVmax = 2n + Z = 8, 10, 12, . . . 18. 
The results are well converged in large model spaces with 
A'max = 18. Let us now compare the ^He calculations in 
the CCSD approximation with the EIHH results. Fig- 
ure [ijb) shows the LITs for F = 10 MeV as a function 
of luq. The CCSD calculation is based on model spaces 
with A'niax — 18 and oscillator frequency hfl = 20 MeV 
and agrees very well with the EIHH result. In these large 
model spaces, the CCSD results are practically indepen- 
dent of the oscillator frequency. We attribute the small 
differences with respect to the exact result to the trun- 
cation inherent in the CCSD approximation. 

It is also interesting to compare dipole response func- 
tions, i.e. inversions of the LIT, obtained from the EIHH 
and the CC methods. For the inversion, we use the min- 
imization procedure of Refs. PB1I12], and require the so- 
lutions of this ill-posed problem to be zero below the 
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about w = 30 and 50 MeV. 
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FIG. 1. (Color online) Convergence of the LIT with F = 
10 MeV as a function of Nmax for an harmonic oscillator 
frequency of TiQ, — 20 MeV (a). Comparison of the LIT at 
F = 10 calculated in the CCSD scheme with A^max = 18 and 
no. = 20 MeV against the LIT from EIHH (b). 



threshold energy wth- The latter is the difference be- 
tween the binding energy of ''He and '^H, being 7 + "^He 
"^H + p the first open reaction channel in the inclu- 
sive process. For the N'^LO tvifo-body potential the bind- 
ing energies are 23.97 (7.37) MeV obtained within the 
CCSD approximation and the particle-removed equation- 
of-motion method [33 , and 25.39 (7.85) MeV with EIHH 
for ''He ('^H), respectively. Figurel2]shows the comparison 
of the response functions obtained by inverting the LIT 
from the CC method and from EIHH. The exact thresh- 
old energy cjth — 17.54 MeV has been used in all the 
inversions. For the CCSD calculations we used a model 
space with iVmax = 18 and hO. = 20 MeV. We found 
that the inversions are insensitive to 7if2. The band for 
the CCSD results (particularly visible at the maximum 
around w ~ 24 MeV) results from inverting the LIT with 
r = 10 MeV and F = 20 MeV, respectively, and esti- 
mates the numerical uncertainty associated with the in- 
version. For the EIHH results, the inversions of the LIT 
with F = 10 and 20 MeV overlap. The CCSD response 
function is close to the virtually exact EIHH result. Ap- 
parently, the small deviations between the CCSD and the 
exact result for the LIT in Fig. translates into small 
deviations in the response function for energies between 
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FIG. 2. (Color online) Comparison of the *He dipole response 
function calculated with CCSD {hQ = 20 MeV and iVmax = 
18) with the EIHH result. 

The ^He benchmark suggests that the CC method can 
be employed for the computation of the dipole response, 
and that theoretical uncertainties with respect to the 
model space and the inversion of the LIT are well con- 
trolled. This makes it exciting to apply the CC method 
for the computation of the giant dipole resonance of ^^O. 

Extension to ^^O. — We employ the CC method for 
the computation of the LIT. Figure [sjja) shows the con- 
vergence of the LIT for F = 10 MeV and a sequence of 
model spaces with hil = 26 MeV and A^max from 8 up to 
18. The convergence is good, and small differences be- 
tween A'lnax = 16 and Nj^ax = 18 are only visible close to 
the maximum. We recall that the convergence of the LIT 
also depends on the choice of the resolution parameter F, 
and that a faster convergence is achieved for increasing 
values of F. In what follows, we show the smallest F for 
which we can obtain a good convergence. 

It is interesting to compare our result of the LIT with 
data. The experiment by Ahrens et al. [S^ measured 
the total photo- absorption cross section a-y on an Oxy- 
gen target with natural abundance (99.762% '^^O) with 
an attenuation method. In the unretarded dipole ap- 
proximation, the cross section is related to the response 
fimction S{uj) of Eq. ([I]) via (a is the fine structure con- 
stant) 



(9) 



We extract S{Ld) from the cross-section data, and apply 
the LIT of Eq. ([s]) to the experimental S{(jj) (which is 
possible because the data cover a broad energy range). 
This allows us to directly compare with the theoretical 
result for the LIT, avoiding the inversion procedure. Fig- 
ure [3jb) shows the comparison between data and the 
CC results obtained in a model space with A'max = 18. 
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FIG. 3. (Color online) Convergence of L(ajo,r) at F = 
10 MeV as a function of Nmax for an harmonic oscillator 
frequency of HQ = 26 MeV (a). Comparison of the LIT at 
F = 10 MeV within CCSD for Nmax = 18 and the Lorentz 
integral transform of Ahrens et al. data W (b). 

The band spanned by the experimental data results from 
Lorentz-transforming the data with the corresponding 
error bars. The band of the CCSD results is obtained 
by varying the harmonic oscillator frequency between 
hVt = 20 and 26 MeV. Even for F = 10 MeV a reso- 
nant shape of the transformed S{ui) data (see them in 
Fig. |4]) is clearly visible. The theoretical and experimen- 
tal results agree within the uncertainties in almost all 
the Wo range. Only from cjq = to about 15 MeV, the 
theory slightly underpredicts the data. Recall that the 
(normalized) Lorentzian kernel is a representation of the 
(5- function. Therefore the integral in of L{ujq, F) is the 
same as the integral in uj of S{uj). Also peak positions are 
approximately conserved. Thus, from the agreement of 
the LITs we can already anticipate that the theoretical 
results well reproduce the centroid of the experimental 
dipole response and the bremsstrahlung sum rule. 

Finally, we perform the inversion of the computed LIT 
and confront theory directly with the data. For the in- 
version, we make the following ansatz for the response 
function 
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FIG. 4. (Color online) Comparison of the ^^O dipole response 
calculated in the CCSD scheme against experimental data by 
Ahrens et al. |39| (triangles with error bars), and Ishkhanov 
et al. [30] (red circles). 



Here, the exponential prefactor is a Gamow factor, and 
/i ~ ^^m„ is the reduced mass with m„ being the nu- 
cleon mass. The charges Zi =1 and Z2 = 1 correspond 
to the first disintegration channel (proton separation) in 
the 7 + ^^O reaction. The binding energy of ^^O in 
the CCSD approximation is 107.24 MeV and the corre- 
sponding threshold energy computed using the particle- 
removed equation-of-motion theory is 14.25 MeV. Note 
that the inversion is sensitive to numerical noise, and 
that the CCSD results using the Harmonic Oscillator ba- 
sis are (only) converged at a few-percent level, and not at 
an ideal sub-percent level. In order to improve the con- 
vergence and make the inversion more stable we employ a 
basis of bound and discretized continuum states obtained 
from diagonalizing a spherical Woods-Saxon potential in 
a discrete plane-wave basis of 35 mesh points for the pro- 
ton and neutron (^5/27 ■Si/27 ^^3/2 partial waves |41]. For 
the remaining partial waves we employed Nmax = 16 
harmonic oscillator shells with hTl = 26 MeV. The in- 
version determines the coefficients Ci and the non-linear 
parameter /3 of the N basis functions by a least-square fit, 
and the results exhibit a corresponding uncertainty. Fig- 
ure[4]shows the CCSD dipole response function and com- 
pares it to the data of Ahrens et al. [S^ , and also to the 
more recent evaluation by Ishkhanov et al. [30", . The 
theoretical band is obtained by inverting two Lorentzian 
widths F = 10, 20 MeV, and by varying the number N 
of basis functions employed in the inversion from 5 to 9. 
This band is a lower estimate of the theoretical error, as it 
does not account e.g. for missing triples in the cluster ex- 
pansion. It is very exciting that the position of the GDR 
in ^^O is nicely explained by our first principle calcula- 
tion whose only ingredient is a NN interaction which fits 
NN scattering data. We observe that the form of the- 
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oretical result is somewhat smeared compared to data. 
The tail region between 40 and 100 MeV is reproduced 
within uncertainties. When integrating the theoretical 
photo-absorption cross section up to 100 MeV we obtain 
an enhancement k — 0.57 — 0.58 of the Thomas-Reiche- 
Kuhn sum rule [59.74^MeV mb(l + k)] . 

Summary. — We presented an ab-initio calculation of 
the GDR in ^^O. Starting from a NN interaction from 
chiral effective field theory that reproduces NN scatter- 
ing data with high accuracy, we computed the Lorentz 
integral transform of the dipole response function with 
the coupled-cluster method. An inversion of the trans- 
form yields the dipole response function. The theoretical 
result reproduces the experimental centroid as well as 
the total dipole response strength (bremsstrahlung sum 
rule) very well, and exhibits a somewhat smeared shape 
compared to the data. This is a first step toward the 
ah initio description of the electromagnetic responses of 
medium-heavy stable nuclei and rare isotopes. 
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